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ABSTRACT 

Late-type low surface brightness galaxies (LSBs) are faint disk galaxies with central maximum stellar sur- 
face densities below lOOMopc"^. The currently favored scenario for their origin is that LSBs have formed in 
fast-rotating halos with large angular momenta. We present the first numerical evidence for this scenario using 
a suite of self-consistent hydrodynamic simulations of a 2.3 x 10^ ^M© galactic halo, in which we investigate 
the correlations between the disk stellar/gas surface densities and the spin parameter of its host halo. A clear 
anti-correlation between the surface densities and the halo spin parameter A is found. That is, as the halo spin 
parameter increases, the scale radius at which the stellar surface density drops below 0.1 Mopc"^ monoton- 
ically increases, while the average stellar surface density of the disk within that radius decreases. The ratio 
of the average stellar surface density for the case of A = 0.03 to that for the case of A = 0.14 reaches more 
than 15. We demonstrate that the result is robust against variations in the baryon fraction, confirming that the 
angular momentum of the host halo is an important driver for the formation of LSBs. 

Subject headings: galaxies: formation — galaxies: evolution — stars: formation — cosmology: dark matter 



1. INTRODUCTION 

Still shrouded in mystery are the origin and true nature of 
dark matter that has long evaded astronomical observations. 
Astrophysical interest in dark matter has been escalated not 
only because it is inferred to dominate the matter mass bud- 
get of the Universe, but because of the simplicity of its dy- 
namics. In particular, dark matter is thought to not only ex- 
ist inside galactic- scale halos, but even dominate the galactic 
gravitational dynamics. Naturally, the records of dark matter 
dynamics throughout the formation of the galactic halos are 
imprinted in the morphology of the present-day galaxies. We 
may therefore advance our understanding of the gravity and 
the nature of dark matter by investigating the galaxies' struc- 
ture and formation history. For this purpose, late-type low 
surface brightness galaxies (LSBs) provide us with a unique 
and intriguing laboratory. 

Late-type LSBs typically refer to the disk galaxies that are 
fainter than the night sky, exhibiting 5-band central surface 
brightnesses /i^^o > 22.7 mag arcsec"^ (Freeman 1970) or 
equivalently central maximum stellar surface densities o < 
lOOMopc-^ (McGaugh et al. 2001). They are different from 
the under-luminous dwarf spheroidals in that the late-type 
LSBs may have large disks and high intrinsic luminosities in 
spite of their low surface brightnesses. Their abundance has 
long been underestimated, but observations now imply that 
LSBs may represent a significant fraction of the galaxy pop- 
ulation (e.g. Impey & Bothun 1997, and references therein). 
Since it was realized that a wealth of LSBs inhabit both cluster 
and field regions, and they may provide an important clue for 
the galaxy formation process (e.g. Impey et al. 1988; Irwin et 
al. 1990; McGaugh et al. 1995; Sprayberry et al. 1996, 1997), 
many studies have explored their physical properties, dynam- 
ical states, host halo structures, spatial distribution, and envi- 
ronmental dependences (e.g. McGaugh & Bothun 1994; Mo 
et al. 1994; Zwaan et al. 1995; de Blok et al. 1996; Mihos et al. 
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1996; Gerritsen & de Blok 1999; Bell et al. 2000; Matthews 
& Wood 2001; Bergmann et al. 2003; Kuzio de Naray et al. 
2004; Matthews et al. 2005; Boissier et al. 2008; Rosenbaum 
et al. 2009; Gao et al. 2010; Galaz et al. 2011; MorelH et al. 
2012; Zhong et al. 2012; CeccarelH et al. 2012). 

While the origin of LSBs is yet to be fully understood, theo- 
retical modelings to comprehend how LSBs have formed par- 
alleled the above phenomenological efforts. The simplest, but 
currently favored scenario states that the difference of LSBs 
from the high surface brightness galaxies (HSBs) results from 
the fast-spinning motion of their host halos (see §2; Dalcanton 
et al. 1997; Jimenez et al. 1997; Mo et al. 1998). To quanti- 
tatively compare the predictions of this theory with the ob- 
served properties of LSBs, Jimenez et al. (1998) utilized sim- 
plified models for the physical, chemical, and spectrophoto- 
metric evolution of disk galaxies, and found that the theory 
nicely reproduces the observed properties of LSBs. Their re- 
sults were later confirmed by Boissier et al. (2003) who used 
larger and higher quality LSB samples (see §2). 

Because LSBs are more dark matter dominant than HSBs 
(e.g. de Blok & McGaugh 1997; Pickering et al. 1997; Mc- 
Gaugh et al. 2001), LSBs have also been employed to test 
the nature of dark matter and gravity. Studies have claimed 
that the rotation curves of LSBs could be useful in constrain- 
ing the properties of dark matter and even modified Newto- 
nian dynamics (e.g. McGaugh & de Blok 1998a,b; de Blok 
& McGaugh 1998; Spergel & Steinhardt 2000; S waters et al. 
2010; Kuzio de Naray & Spekkens 2011). Recently Lee et 
al. (2012) speculated that the abundance of LSBs could be a 
testbed for modified gravity. Using a high-resolution simula- 
tion, they demonstrated that modified gravity may spin up the 
low mass galactic halos and significantly alter the abundance 
of LSBs. Nevertheless, in order to utilize LSBs as a probe 
into the nature of dark matter and gravity, it is indispensable 
to first properly constrain the driver that produces LSBs. Al- 
though the galactic evolution models of Jimenez et al. (1998) 
and Boissier et al. (2003) supported the hypothesis that the 
LSBs form in the halos with large angular momenta, their re- 
sults were subject to the simplified prescriptions for the for- 
mation of disk galaxies. To improve the situation and advance 



2 



J. Kim & J. Lee 




^0 100 200 300 400 100 200 300 400 

X (70 kpc) X (70 kpc) 

Fig. 1. — Face-on stellar surface densities in a central 70 kpc box at 1.49 Gyrs after the start of the simulation for four different halo spin parameters, 
A = 0.03, 0.06, 0.10, 0.14 in the top-left, top-right, bottom-left, and bottom-right panel, respectively. The initial baryonic mass fraction in the halo is /b = 0.10. 
Stars of age less than 1 Gyr are used to estimate surface densities. More information on the suite of simulations and analysis is provided in §3 and §4.1. 



our understanding of the origin of LSBs, a comprehensive, 
self-consistent hydrodynamic simulation is greatly needed. 

In light of this development, we for the first time numeri- 
cally examine whether or not the angular momentum of the 
host halo is an important driver for determining the galactic 
surface brightness. Using a suite of comprehensive hydrody- 
namic simulations, we quantitatively investigate the correla- 
tions between the surface density of a disk galaxy and the spin 
parameter of its host halo without having to rely on a simpli- 
fied prescription for galactic evolution. The outline of this 
article is as follows. We briefly describe the currently popu- 
lar model for the origin of LSBs in §2, and the set-up of our 
numerical experiments in §3. The analysis of our simulation 
results is presented in §4 focusing on the correlation between 
the stellar/gas surface densities and the halo spin parameter. 
The robustness of our result is also examined against varia- 
tions in the initial baryonic mass fraction. Finally we discuss 
the implications of our results and draw a conclusion in §5. 

2. A BRIEF REVIEW OF THE MODEL FOR THE ORIGIN OF LOW 
SURFACE BRIGHTNESS GALAXIES 

A dark matter halo acquires its angular momentum J 
through the tidal interactions with the surrounding matter dis- 
tribution at its proto-halo stage (Peebles 1969; Doroshkevich 



1970; White 1984). Since the magnitude of the halo an- 
gular momentum is dependent upon its virial mass Myir, it 
is often convenient to adopt a dimensionless spin parame- 
ter A = |J|/(2GMvir^vir)^^^ with the halo virial radius 7?vir to 
quantify how fast a halo rotates around its axis of symmetry. 
Simulations have shown that the probability density of X is 
approximated by a log-normal distribution in ACDM cosmol- 
ogy, being relatively insensitive to mass scale, environment, 
and redshift (e.g. Bullock et al. 2001). 

In the gravitationally self-consistent model for the forma- 
tion of disk galaxies, a rotationally- supported disk forms in a 
potential well of a dark matter halo through a dissipative grav- 
itational collapse of the baryonic content (Fall & Efstathiou 
1980). The resulting baryonic disk shares the same tidally- 
induced specific angular momentum j (angular momentum 
per unit mass) with its host halo. Based on this model, Dal- 
canton et al. (1997) proposed a scenario that the large specific 
angular momentum of a dark matter halo disperses the gas to 
a wider extent, leading to a large disk size and thus a low gas 
surface density (see also Jimenez et al. 1997; Mo et al. 1998). 
To relate the gas surface density to the stellar surface density 
and to perform a systematic test of this scenario against ob- 
servations, Jimenez et al. (1998) computed various galactic 
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Fig. 2. — Same as Figure 1 but from the edge-on angle. A = 0.03, 0.06, 0.10, 0.14 in the top-left, top-right, bottom-left, and bottom-right panel, respectively. 



properties including surface brightnesses by adopting the fol- 
lowing prescriptions for galaxy evolution: (a) 3.n isothermal 
spherical halo, (b) 3. model for gas infall rate that reproduces 
the observed properties of the Milky Way, fcj the Schmidt law 
for star formation (Schmidt 1959), and (d) the galactic chem- 
ical evolution model by Matteucci & Francois (1989). They 
found nice agreement between the resulting predictions and 
the observed properties of LSBs such as surface brightness, 
metallicity, and colors, provided that the host halo satisfies the 
condition of A > 0.06. This work was extended by Boissier 
et al. (2003) to a larger set of observational data and a wider 
range of halo spin parameters. They drew the same conclu- 
sion that a LSB is likely to form in a dark matter halo with a 
large spin parameter. They also noted that some modulation 
in the star formation rate history would be desired to achieve 
a better agreement between the theory and observations. 

However, the conclusion reached by Jimenez et al. (1998) 
and Boissier et al. (2003) was contingent upon how the gas in- 
fall rate and star formation rate were prescribed in their overly 
simplified evolution models. It is still inconclusive whether 
the theoretical model by Dalcanton et al. (1997) would hold 
robustly in reality where the gas infall and star formation can 
hardly be specified by a simple analytic recipe. In what fol- 
lows, we describe a suite of comprehensive, first-principles 



numerical simulations of a galactic halo with realistic initial 
conditions, designed to better address this long-standing prob- 
lem. Most importantly, unlike previous studies we do not rely 
on simplified prescriptions for star formation and gas infall 
rates, which for the first time enables us to examine the ro- 
bustness of the theory in the most realistic set-up yet. 

3. INITIAL CONDITIONS AND PHYSICS IN THE CODE 

We simulate the evolution of a dwarf- sized galactic halo 
with varying halo spin parameters and initial baryon fractions. 
The initial conditions of our experiment and the physics in- 
cluded in the code are explained in this section. 

We set up an isolated halo of total mass 2.3 x 10^^ 
with four halo spin parameters, A = 0.03,0.06,0.10,0.14, 
and with two different initial baryon fractions for each A, 
/b = (M^ +Mgas)/Mtotai = 0.05 and 0.10 (see Table 1; 95% 
dark matter + 5% gas for /b =0.05, and 90% dark matter + 
10% gas for /b = 0.10). We first create a dataset of 10^ col- 
lisionless particles, to which we add gas by splitting the par- 
ticles (with an initial metalliticy of 0.003 Z©). Gas and dark 
matter follow the same shapes of a Navarro-Frenk- White pro- 
file (Navarro et al. 1997) with concentration c = 10, and have 
the same value of A as an averaged circular motion. This way, 
as the gas cools down, the progenitor forms a disk galaxy em- 
bedded in a gaseous halo. Note that this procedure is very dif- 
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Fig. 3. — Azimuthally-averaged stellar surface densities at 1.49 Gyrs after 
the start of the simulation with /b = 0.10 for A = 0.03, 0.06, 0.10, 0.14 as a 
dashed, solid, dot-dashed, and dotted line, respectively. The scale radius 
at which the stellar surface density falls below 0.1 Mq pc~^ monotonically 
increases as the halo spin parameter increases (see §4.1 and Table 1). 



ferent from that of Jimenez et al. (1998) where a galactic disk 
already exists when their calculation starts. For detailed de- 
scriptions to set up the initial condition, see Kim et al. (2012). 

We then follow the evolution of each galactic halo in a 1 
Mpc^ box for ^ 2 Gyrs using the ZEUS hydrodynamics solver 
included in the publicly available adaptive mesh refinement 
Enzo-2.1? Multi- species chemistry (H, He, He^, He^^, 
e~) and non-equilibrium cooling modules are employed to 
calculate radiative losses. Cooling by metals is estimated in 
gas above 10^ K with the rates tabulated by Sutherland & Do- 
pita (1993), and below 10"^ K by Glover & Jappsen (2007). 
We refine the cells by factors of two in each axis on the over- 
densities of gas and dark matter. In the finest cell of size 
Ax = 122 pc, a star cluster particle is produced with initial 
mass M^^^ = O.SpgasAx^ when (a) the proton number density 
exceeds 3.2 cm~^, (b)thQ velocity flow is converging, fcj the 
cooling time ^cooi is shorter than the dynamical time ^dyn of 
the cell, and (d) the particle produced has at least 104000 Mq. 
For each star cluster particle, 1.5 x 10~^ of its rest mass en- 
ergy and 20% of its mass are returned to the gas phase over 12 
^dyn- This represents various types of stellar feedback, chiefly 
the injection of thermal energy by supernova explosion. More 
information on adopted baryonic physics can be found in Kim 
et al. (2009) and Kim et al. (201 1). 

4. SIMULATION RESULTS 

We now describe the results of our simulations with varying 
halo spin parameters and baryon fractions. First, we focus on 
the correlations between surface densities of a galactic disk 
and the spin parameter of its host halo in the simulations with 
an initial baryonic mass fraction of 10%. 

4.1. Galactic Disks with Varying Spin Parameters 

We start by calculating the stellar surface density of the 
simulated disk galaxy for each case of A at 1.49 Gyrs af- 
ter the start of the simulation with an initial baryon fraction 
/b = 0.10, the ratio of baryonic mass to total mass. Shown 
in Figures 1 and 2 are the snapshots of face-on and edge-on 

^ http://enzo-project.org/ 
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Fig. 4. — Same as Figure 4 but for azimuthally-averaged gas surface den- 
sities. The size of a gas disk increases with the halo spin parameter. 



Stellar surface densities using the stars of age less than 1 Gyr. 
In each figure, the four panels correspond to the four differ- 
ent values of halo spin parameters, A = 0.03, 0.06, 0.10, 0.14. 
The snapshots are generated in a 70 kpc box with uniform 150 
pc resolution which can be regarded as an equivalent of an 
aperture size in photometric observations. We also point out 
that, in these figures and the subsequent analysis, we exclude 
the stars formed in the first 0.49 Gyrs since they would most 
possibly reflect the artificial burst of star formation in the ini- 
tially unimpeded collapse, caused by an unstable gas profile 
at the beginning of the simulation.^ From these figures, it is 
clearly noticeable that the size of the stellar disk is larger in 
the case of a higher value of A, indicating that the large an- 
gular momentum of a dark matter halo disperses the baryonic 
content to a wider extent. 

To quantify the change in the disk size resulting from dif- 
ferent halo spin parameters, the radial profile of the stellar 
surface density, I^^(r) (in the unit of MqPC~^), is determined 
by taking the azimuthal average of the stellar density on the 
disk plane at 1.49 Gyrs after the start of the simulation. The 
result of this exercise is displayed in Figure 3. Similarly, the 
radial profile of the azimuthally-averaged gas surface den- 
sity, Lgas(r), is shown in Figure 4. Then, the scale radius 
Ts is defined as the radius at which the stellar surface den- 
sity drops below 0.1 Mopc"^. Table 1 lists the scale radii for 
the four different A , and the average stellar surface densities 
within those radii, Z^,ave- One can observe that the scale ra- 
dius Ts monotonically increases as A increases, from 8.1 kpc 
at A = 0.03 to 32.5 kpc at A = 0. 14. The net effects are an en- 
larged photometric size of the stellar disk, and a lowered aver- 
age stellar surface density Z^,ave- 278.2 Mopc"^ at A = 0.03 
to 15.2M0pc~^ at A = 0.14. As a result, the ratio of E^,ave 
for the case of A = 0.03 to that of A = 0.14 becomes larger 
than 15. This result confirms the idea that the faster rota- 
tion of the host halo disperses the disk gas to a larger extent, 
which then leads to a lower stellar surface density averaged on 
the disk. Assuming a plausible stellar mass-to-light ratio, the 
low stellar surface density is translated into the low surface 

^ It should also be noted that we arbitrarily choose the epoch of 1 .49 Gyrs 
for the analysis. Excluding the stars formed in the first 0.49 Gyrs we, there- 
fore, implicitly assume that the stars younger than 1 Gyr form the disk of a 
dwarf- sized galaxy. While 1 Gyr may correspond to an interval between ma- 
jor disturbances for dwarf-sized galaxies, our choice is certainly very fiducial. 
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Fig. 5. — Radial profiles of enclosed baryonic mass (stars+gas) at 1.49 Gyrs 
after the start of the simulation with /b = 0.10 for A = 0.03, 0.06, 0.10, 0.14 
as a thick dashed, solid, dot-dashed, and dotted line, respectively. The initial 
radial profile is also shown for comparison as a thin dotted line. While all the 
halos show the results of dissipative collapses towards the center, the halo of 
A = 0.03 has the densest core in the baryon distribution. 

brightness, giving consistent results with the previous studies 
by Jimenez et al. (1998) and Boissier et al. (2003). 

It is also very informative to investigate how the mass pro- 
files are correlated with the halo spin parameter. We evaluate 
the profiles of baryonic mass enclosed within a sphere of ra- 
dius r from the galactic center of mass at 1 .49 Gyrs after the 
start of the simulation. The resulting radial profiles in Fig- 
ure 5 allow us to cleanly see the discrepancy across different 
runs in a cumulative manner. While all the halos show the 
results of dissipative collapses towards the center, the halo of 
the lowest spin parameter A = 0.03 has the densest core in the 
baryonic distribution, enclosing most of the baryons within 
~ 3 kpc from the center. We also inspect the half-mass ra- 
dius that encloses a half of the total baryonic mass inside the 
halo, ~ 2.2 X lO^^M©. As is anticipated, the half-mass ra- 
dius increases almost monotonically with A, from 1.4 kpc at 
A = 0.03 to 8.0 kpc at A = 0.14. Finally Figure 6 demon- 
strates the profiles of the dark matter mass enclosed within a 
sphere of radius r. One can easily notice that the larger the 
halo spin parameter is, the less dense the dark matter distri- 
bution becomes in the core region (r < 1 kpc). In particu- 
lar, the mass enclosed within a 1 kpc sphere in the A =0.03 
case is ~ 2.8 X lO^M©, whereas that of A = 0.14 case is only 
~ 6.0 X lO^M©. This implies that dark matter itself is also 
gradually dispersed outwards by the fast rotation of the halo. 

4.2. Robustness Against Variations in the Baryon Fraction 

Now that we have found a strong correlation between the 
halo spin parameter and the stellar/gas surface densities of a 
galactic disk, it will be very illuminating to examine whether 
or not this correlation is robust against variations in the initial 
baryonic mass fraction /b. This examination may also allow 
us to determine which factor is more dominant in determining 
the scale radius and galactic surface density, between the ab- 
solute amount of gas and the angular momentum of the host 
halo. Recall that the results presented in §4.1 are all obtained 
by the simulations with /b = 0.10. 

We here repeat the same analysis in Figure 3 but on a 
suite of simulations with a different initial baryon fraction, 
/b = 0.05, and plot the result in Figure 7. First, as is ob- 
viously expected, a decreased value of /b also reduces the 
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Fig. 6. — Same as Figure 5 but of the enclosed dark matter mass. Note that 
the halo of A = 0.03 exhibits the most centrally-concentrated density profile. 
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Fig. 7. — Same as Figure 3 but for the case of the initial baryonic mass 
fraction /b = 0.05 instead of /b = 0.10. The trend of an increasing scale 
radius with a larger halo spin parameter persists. 



overall stellar surface density for a halo with the same spin 
parameter. What is more intriguing is that the positive corre- 
lation between the scale radius r^ (see §4.1) and the halo spin 
parameter A still holds strong in simulations with a lower /bi 
a faster rotation of the host halo induces a larger scale radius 
of the stellar disk. Table 1 lists the scale radii and the average 
stellar surface densities of the galactic disks within those radii 
for /b =0.05, which cleanly shows that the correlations be- 
tween Ts and A holds the same for both /b's. Interestingly, the 
difference between two rs's for /b = 0.05 and 0.10 with the 
same A is surprisingly small. For example, the scale radius 
of a disk in the halo of A = 0.10 and /b = 0.10 is 26.8 kpc, 
only slightly expanded from 26.0 kpc in the halo of the same 
A but of /b = 0.05. This is despite the fact that the initial gas 
supply is boosted by a factor of two. Although it should be 
considered provisional this result may suggest that, between 
A and /b, the more influential factor for determining the scale 
radius is A, the spin parameter of the halo. 

4.3. Notes on the Performed Runs 

Before we proceed to conclude, it is worth to note a few 
points on the simulations we have described so far. We draw 
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TABLE 1 

HALO SPIN PARAMETER, DISK SCALE RADIUS, AND AVERAGE STELLAR SURFACE DENSITY^ 





/b = 


0.05 




/b = 


0.10 


halo spin parameter A 


scale radius [kpc] 


^★,ave 


[Mopc-2] 


scale radius [kpc] 


^*,ave [Mopc-^] 


0.03 


7.2 




157.0 


8.1 


278.2 


0.06 


13.3 




49.1 


17.3 


60.2 


0.10 


26.0 




11.7 


26.8 


24.6 


0.14 


32.2 




6.5 


32.5 


15.2 



^For definitions and detailed explanation, see §4.1. 

the reader's attention to the central maximum stellar surface 
densities of the simulated galaxies in Figure 3, which still 
remain above or very close to lOOMopc"^ regardless of A. 
Therefore, at face value, even the largest halo spin parameter 
seemingly fails to produce a LSB galaxy in its conventional 
definition (see §1). However, these simulations are under 
a particular assumption of a relatively high baryon fraction 
of /b = 0.10. By comparing Figures 3 and 7, one observes 
that the central peak stellar surface density indeed drops sig- 
nificantly when the initial gas supply is reduced by a fac- 
tor of two. Had the galaxy been simulated in a fast-rotating 
halo (e.g. A > 0.10) with an even smaller baryon fraction 
(e.g. /b < 0.02), it would likely have developed into a LSB 
galaxy with a central maximum stellar surface density less 
than lOOM0pc~^. Further, by running simulations in an iso- 
lated set-up, we also have made an implicit assumption that 
the halo does not experience any interaction with its surround- 
ing environment, such as mergers and harassment. Including 
such effects might have drastically changed the surface densi- 
ties of the simulated galaxies. 

In summary, the reported simulations are aimed to make a 
comparison between the galaxies that differ only by their halo 
spin parameters. Our experiment is however not specifically 
designed to turn any galaxy into a LSB with a large halo spin 
parameter. In reality, the galactic surface brightness is most 
likely determined by the combination of a number of factors 
such as halo angular momentum, baryon fraction, and envi- 
ronmental effects (see §5). 

5. DISCUSSION AND CONCLUSION 

We have presented the first numerical evidence for the hy- 
pothesis that late-type LSBs form in the halos with large an- 
gular momenta. A suite of comprehensive hydrodynamic sim- 
ulations of a 2.3 X 10^^ Mq galactic halo is employed with 
different halo spin parameters and baryon fractions. We have 
investigated the correlations between the stellar/gas surface 



densities of a galactic disk and the spin parameter of its host 
halo. A clear signal of anti-correlation is found between the 
surface densities and the halo spin parameter. That is, as 
the halo spin parameter increases, the scale radius, defined 
as the radius at which the stellar surface density drops below 
0.1 M0pc~^, monotonically increases, while the average stel- 
lar surface density of the disk within that radius decreases. 
The ratio of the average stellar surface density for the case 
of A = 0.03 to that for the case of A = 0.14 is larger than 
15. Even dark matter itself is dispersed outwards by a large 
halo angular momentum. We also have demonstrated that the 
result is robust against variations in the initial value of the 
baryonic mass fraction. 

Our investigation confirms that one of the most impor- 
tant drivers for the formation of LSBs is the spin parame- 
ters of their host halos. This work provides a critical link be- 
tween observations and the theoretical models for the origin 
of LSBs, which has long been missing in the previous studies. 
It also enables us to better comprehend the formation process 
of LSBs, and to utilize LSBs for the purpose of constraining 
the nature of invisible component that dominates galactic dy- 
namics. Even so, the reported calculation alone cannot decide 
whether the halo angular momentum is the most determini- 
tive factor in the formation of LSBs. There is a large multi- 
dimensional parameter space to explore, including the halo 
angular momentum, halo density profile, baryon fraction, gas 
infall rate, and stellar feedback, to name a few. 
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by the Korea government (MEST, No.20 12-00049 16) and 
from the National Research Foundation of Korea to the Center 
for Galaxy Evolution Research. 
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